Three near-complete genome assemblies reveal substantial centromere dynamics from diploid to tetraploid in Brachypodium genus.

BACKGROUND: Centromeres are critical for maintaining genomic stability in eukaryotes, and their turnover shapes genome architectures and drives karyotype evolution. However, the co-evolution of centromeres from different species in allopolyploids over millions of years remains largely unknown. RESULTS: Here, we generate three near-complete genome assemblies, a tetraploid Brachypodium hybridum and its two diploid ancestors, Brachypodium distachyon and Brachypodium stacei. We detect high degrees of sequence, structural, and epigenetic variations of centromeres at base-pair resolution between closely related Brachypodium genomes, indicating the appearance and accumulation of species-specific centromere repeats from a common origin during evolution. We also find that centromere homogenization is accompanied by local satellite repeats bursting and retrotransposon purging, and the frequency of retrotransposon invasions drives the degree of interspecies centromere diversification. We further investigate the dynamics of centromeres during alloploidization process, and find that dramatic genetics and epigenetics architecture variations are associated with the turnover of centromeres between homologous chromosomal pairs from diploid to tetraploid. Additionally, our pangenomes analysis reveals the ongoing variations of satellite repeats and stable evolutionary homeostasis within centromeres among individuals of each Brachypodium genome with different polyploidy levels. CONCLUSIONS: Our results provide unprecedented information on the genomic, epigenomic, and functional diversity of highly repetitive DNA between closely related species and their allopolyploid genomes at both coarse and fine scale.

Group VII ethylene response factors forming distinct regulatory loops mediate submergence responses.

Group VII ethylene response factors (ERFVIIs), whose stability is oxygen concentration-dependent, play key roles in regulating hypoxia response genes in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) during submergence. To understand the evolution of flooding tolerance in cereal crops, we evaluated whether Brachypodium distachyon ERFVII genes (BdERFVIIs) are related to submergence tolerance. We found that three BdERFVIIs, BdERF108, BdERF018, and BdERF961, form a feedback regulatory loop to mediate downstream responses. BdERF108 and BdERF018 activated the expression of BdERF961 and PHYTOGLOBIN 1 (PGB1), which promoted nitric oxide turnover and preserved ERFVII protein stability. The activation of PGB1 was subsequently counteracted by increased BdERF961 accumulation through negative feedback regulation. Interestingly, we found that OsERF67, the orthologue of BdERF961 in rice, activated PHYTOGLOBIN (OsHB2) expression and formed distinct regulatory loops during submergence. Overall, the divergent regulatory mechanisms exhibited by orthologs collectively offer perspectives for the development of submergence-tolerant crops.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum.

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium.

The "genomic shock" hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession "Bhyb26." We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species' similarity in transposable element load may account for the subtlety of the observed genome dominance.

Genome Wide Identification and Annotation of NGATHA Transcription Factor Family in Crop Plants.

The NGATHA (NGA) transcription factor (TF) belongs to the ABI3/VP1 (RAV) transcriptional subfamily, a subgroup of the B3 superfamily, which is relatively well-studied in Arabidopsis. However, limited data are available on the contributions of NGA TF in other plant species. In this study, 207 NGA gene family members were identified from a genome-wide search against Arabidopsis thaliana in the genome data of 18 dicots and seven monocots. The phylogenetic and sequence alignment analyses divided NGA genes into different clusters and revealed that the numbers of genes varied depending on the species. The phylogeny was followed by the characterization of the Solanaceae (tomato, potato, capsicum, tobacco) and Poaceae (Brachypodium distachyon, Oryza sativa L. japonica, and Sorghum bicolor) family members in comparison with A. thaliana. The gene and protein structures revealed a similar pattern for NGA and NGA-like sequences, suggesting that both are conserved during evolution. Promoter cis-element analysis showed that phytohormones such as abscisic acid, auxin, and gibberellins play a crucial role in regulating the NGA gene family. Gene ontology analysis revealed that the NGA gene family participates in diverse biological processes such as flower development, leaf morphogenesis, and the regulation of transcription. The gene duplication analysis indicates that most of the genes are evolved due to segmental duplications and have undergone purifying selection pressure. Finally, the gene expression analysis implicated that the NGA genes are abundantly expressed in lateral organs and flowers. This analysis has presented a detailed and comprehensive study of the NGA gene family, providing basic knowledge of the gene, protein structure, function, and evolution. These results will lay the foundation for further understanding of the role of the NGA gene family in various plant developmental processes.

Brachypodium: 20 years as a grass biology model system; the way forward?

It has been 20 years since Brachypodium distachyon was suggested as a model grass species, but ongoing research now encompasses the entire genus. Extensive Brachypodium genome sequencing programmes have provided resources to explore the determinants and drivers of population diversity. This has been accompanied by cytomolecular studies to make Brachypodium a platform to investigate speciation, polyploidisation, perenniality, and various aspects of chromosome and interphase nucleus organisation. The value of Brachypodium as a functional genomic platform has been underscored by the identification of key genes for development, biotic and abiotic stress, and cell wall structure and function. While Brachypodium is relevant to the biofuel industry, its impact goes far beyond that as an intriguing model to study climate change and combinatorial stress.

MiR398-regulated antioxidants contribute to Bamboo mosaic virus accumulation and symptom manifestation.

Virus infections that cause mosaic or mottling in leaves commonly also induce increased levels of reactive oxygen species (ROS). However, how ROS contributes to symptoms is less well documented. Bamboo mosaic virus (BaMV) causes chlorotic mosaic symptoms in both Brachypodium distachyon and Nicotiana benthamiana. The BaMV △CPN35 mutant with an N-terminal deletion of its coat protein gene exhibits asymptomatic infection independently of virus titer. Histochemical staining of ROS in mock-, BaMV-, and BaMV△CPN35-infected leaves revealed that hydrogen peroxide (H2O2) accumulated solely in BaMV-induced chlorotic spots. Moreover, exogenous H2O2 treatment enhanced yellowish chlorosis in BaMV-infected leaves. Both BaMV and BaMV△CPN35 infection could induce the expression of Cu/Zu superoxide dismutase (CSD) antioxidants at messenger RNA and protein level. However, BaMV triggered the abundant accumulation of full-length NbCSD2 preprotein (prNbCSD2, without transit peptide cleavage), whereas BaMV△CPN35 induced a truncated prNbCSD2. Confocal microscopy showed that majority of NbCSD2-green fluorescent protein (GFP) predominantly localized in the cytosol upon BaMV infection, but BaMV△CPN35 infection tended to cause NbCSD2-GFP to remain in chloroplasts. By 5'-RNA ligase-mediated rapid amplification of cDNA ends, we validated CSDs are the targets of miR398 in vivo. Furthermore, BaMV infection increased the level of miR398, while the level of BaMV titer was regulated positively by miR398 but negatively by CSD2. In contrast, overexpression of cytosolic form NbCSD2, impairing the transport into chloroplasts, greatly enhanced BaMV accumulation. Taken together, our results indicate that induction of miR398 by BaMV infection may facilitate viral titer accumulation, and cytosolic prNbCSD2 induction may contribute to H2O2 accumulation, resulting in the development of BaMV chlorotic symptoms in plants.

Genes involved in mRNA surveillance are induced in Brachypodium distachyon under cadmium toxicity.

BACKGROUND: Cd accumulation in plant cells results in dramatic problems including oxidative stress and inhibition of vital enzymes. It also affects mineral uptakes by disrupting membrane permeability. Interaction among Cd and other plant nutrient elements changes the nutritional contents of crops and reduces their yield. METHODS AND RESULTS: In the present study, Cd stress in Brachypodium distachyon led to the upregulation of some heavy metal transport genes (influx or efflux) encoding cation-efflux proteins, heavy metal-associated proteins and NRAMP proteins. The Arabidopsis orthologs of the differentially expressed B. distachyon genes (DEGs) under Cd toxicity were identified, which exhibited Bradi4g26905 was an ortholog of AtALY1-2. Detailed co-expression network and gene ontology analyses found the potential involvement of the mRNA surveillance pathway in Cd tolerance in B. distachyon. These genes were shown to be downregulated by sulfur (S) deficiency. CONCLUSIONS: This is the first transcriptomic study investigating the effect of Cd toxicity in B. distachyon, a model plant for genomic studies in Poaceae (Gramineae) species. The results are expected to provide valuable information for more comprehensive research related to heavy metal toxicity in plants.

In Search of Monocot Phosphodiesterases: Identification of a Calmodulin Stimulated Phosphodiesterase from Brachypodium distachyon.

In plants, rapid and reversible biological responses to environmental cues may require complex cellular reprograming. This is enabled by signaling molecules such as the cyclic nucleotide monophosphates (cNMPs) cAMP and cGMP, as well as Ca2+. While the roles and synthesis of cAMP and cGMP in plants are increasingly well-characterized, the "off signal" afforded by cNMP-degrading enzymes, the phosphodiesterases (PDEs), is, however, poorly understood, particularly so in monocots. Here, we identified a candidate PDE from the monocot Brachypodium distachyon (BDPDE1) and showed that it can hydrolyze cNMPs to 5'NMPs but with a preference for cAMP over cGMP in vitro. Notably, the PDE activity was significantly enhanced by Ca2+ only in the presence of calmodulin (CaM), which interacts with BDPDE1, most likely at a predicted CaM-binding site. Finally, based on our biochemical, mutagenesis and structural analyses, we constructed a comprehensive amino acid consensus sequence extracted from the catalytic centers of annotated and/or experimentally validated PDEs across species to enable a broad application of this search motif for the identification of similar active sites in eukaryotes and prokaryotes.

Programmed Cell Death in Developing Brachypodium distachyon Grain.

The normal developmental sequence in a grass grain entails the death of several maternal and filial tissues in a genetically regulated process termed programmed cell death (PCD). The progression and molecular aspects of PCD in developing grains have been reported for domesticated species such as barley, rice, maize and wheat. Here, we report a detailed investigation of PCD in the developing grain of the wild model species Brachypodium distachyon. We detected PCD in developing Brachypodium grains using molecular and histological approaches. We also identified in Brachypodium the orthologs of protease genes known to contribute to grain PCD and surveyed their expression. We found that, similar to cereals, PCD in the Brachypodium nucellus occurs in a centrifugal pattern following anthesis. However, compared to cereals, the rate of post-mortem clearance in the Brachypodium nucellus is slower. However, compared to wheat and barley, mesocarp PCD in Brachypodium proceeds more rapidly in lateral cells. Remarkably, Brachypodium mesocarp PCD is not coordinated with endosperm development. Phylogenetic analysis suggests that barley and wheat possess more vacuolar processing enzymes that drive nucellar PCD compared to Brachypodium and rice. Our expression analysis highlighted putative grain-specific PCD proteases in Brachypodium. Combined with existing knowledge on grain PCD, our study suggests that the rate of nucellar PCD moderates grain size and that the pattern of mesocarp PCD influences grain shape.

BdFAR4, a root-specific fatty acyl-coenzyme A reductase, is involved in fatty alcohol synthesis of root suberin polyester in Brachypodium distachyon.

Suberin is a complex hydrophobic polymer of aliphatic and phenolic compounds which controls the movement of gases, water, and solutes and protects plants from environmental stresses and pathogenic infection. The synthesis and regulatory pathways of suberin remain unknown in Brachypodium distachyon. Here we describe the identification of a B. distachyon gene, BdFAR4, encoding a fatty acyl-coenzyme A reductase (FAR) by a reverse genetic approach, and investigate the molecular relevance of BdFAR4 in the root suberin synthesis of B. distachyon. BdFAR4 is specifically expressed throughout root development. Heterologous expression of BdFAR4 in yeast (Saccharomyces cerevisiae) afforded the production of C20:0 and C22:0 fatty alcohols. The loss-of-function knockout of BdFAR4 by CRISPR/Cas9-mediated gene editing significantly reduced the content of C20:0 and C22:0 fatty alcohols associated with root suberin. In contrast, overexpression of BdFAR4 in B. distachyon and tomato (Solanum lycopersicum) resulted in the accumulation of root suberin-associated C20:0 and C22:0 fatty alcohols, suggesting that BdFAR4 preferentially accepts C20:0 and C22:0 fatty acyl-CoAs as substrates. The BdFAR4 protein was localized to the endoplasmic reticulum in Arabidopsis thaliana protoplasts and Nicotiana benthamiana leaf epidermal cells. BdFAR4 transcript levels can be increased by abiotic stresses and abscisic acid treatment. Furthermore, yeast one-hybrid, dual-luciferase activity, and electrophoretic mobility shift assays indicated that the R2R3-MYB transcription factor BdMYB41 directly binds to the promoter of BdFAR4. Taken together, these results imply that BdFAR4 is essential for the production of root suberin-associated fatty alcohols, especially under stress conditions, and that its activity is transcriptionally regulated by the BdMYB41 transcription factor.

A novel wheat ASR gene, TaASR2D, enhances drought tolerance in Brachypodium distachyon.

Abscisic acid-, stress-, and ripening-induced (ASR) proteins play an important role in protecting plants against adverse environmental conditions. Here, we identified 24 ASR genes in the wheat genome and analyzed their characteristics. Among these, five ASR genes highly induced by abscisic acid (ABA) and polyethylene glycol were cloned and further characterized. The TaASR genes were expressed in response to different abiotic stresses and ABA and were found to be localized in the nucleus and plasma membrane of transformed tobacco cells. Brachypodium distachyon transgenic plants overexpressing TaASR2D showed enhanced drought tolerance by regulating leaf transpiration. The expression levels of stress-related and ABA-responsive genes were higher in transgenic plants than in wild-type plants under drought stress conditions. Moreover, overexpression of TaASR2D increased the levels of both endogenous ABA and hydrogen peroxide in response to drought stress, and these plants showed hypersensitivity to exogenous ABA at the germination stage. Furthermore, plants overexpressing TaASR2D showed increased stomatal closure. Further analysis revealed that TaASR2D interacts with ABA biosynthesis and stress-related proteins in yeast and tobacco plants. Collectively, these findings indicate that TaASR2D plays an important role in the response of plants to drought stress by regulating the ABA biosynthesis pathway and redox homeostasis system.

The Role of Brachypodium distachyon Wall-Associated Kinases (WAKs) in Cell Expansion and Stress Responses.

The plant cell wall plays a critical role in signaling responses to environmental and developmental cues, acting as both the sensing interface and regulator of plant cell integrity. Wall-associated kinases (WAKs) are plant receptor-like kinases located at the wall-plasma membrane-cytoplasmic interface and implicated in cell wall integrity sensing. WAKs in Arabidopsis thaliana have been shown to bind pectins in different forms under various conditions, such as oligogalacturonides (OG)s in stress response, and native pectin during cell expansion. The mechanism(s) WAKs use for sensing in grasses, which contain relatively low amounts of pectin, remains unclear. WAK genes from the model monocot plant, Brachypodium distachyon were identified. Expression profiling during early seedling development and in response to sodium salicylate and salt treatment was undertaken to identify WAKs involved in cell expansion and response to external stimuli. The BdWAK2 gene displayed increased expression during cell expansion and stress response, in addition to playing a potential role in the hypersensitive response. In vitro binding assays with various forms of commercial polysaccharides (pectins, xylans, and mixed-linkage glucans) and wall-extracted fractions (pectic/hemicellulosic/cellulosic) from both Arabidopsis and Brachypodium leaf tissues provided new insights into the binding properties of BdWAK2 and other candidate BdWAKs in grasses. The BdWAKs displayed a specificity for the acidic pectins with similar binding characteristics to the AtWAKs.

A triticale tapetal non-specific lipid transfer protein (nsLTP) is translocated to the pollen cell wall.

KEY MESSAGE: A Triticeae type III non-specific lipid transfer protein (nsLTP) was shown for the first time to be translocated from the anther tapetum to the pollen cell wall. Two anther-expressed non-specific lipid transfer proteins (nsLTPs) were identified in triticale (× Triticosecale Wittmack). LTPc3a and LTPc3b contain a putative signal peptide sequence and eight cysteine residues in a C-Xn-C-Xn-CC-Xn-CXC-Xn-C-Xn-C pattern. These proteins belong to the type III class of nsLTPs which are expressed exclusively in the inflorescence of angiosperms. The level of LTPc3 transcript in the anther was highest at the tetrad and uninucleate microspore stages, and absent in mature pollen. In situ hybridization showed that LTPc3 was expressed in the tapetal layer of the developing triticale anther. The expression of the LTPc3 protein peaked at the uninucleate microspore stage, but was also found to be associated with the mature pollen. Accordingly, an LTPc3a::GFP translational fusion expressed in transgenic Brachypodium distachyon first showed activity in the tapetum, then in the anther locule, and later on the mature pollen grain. Altogether, these results represent the first detailed characterization of a Triticeae anther-expressed type III nsLTP with possible roles in pollen cell wall formation.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors.

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum.

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Insertional mutagenesis of Brachypodium distachyon using the Tnt1 retrotransposable element.

Brachypodium distachyon is an annual C3 grass used as a monocot model system in functional genomics research. Insertional mutagenesis is a powerful tool for both forward and reverse genetics studies. In this study, we explored the possibility of using the tobacco retrotransposon Tnt1 to create a transposon-based insertion mutant population in B. distachyon. We developed transgenic B. distachyon plants expressing Tnt1 (R0) and in the subsequent regenerants (R1) we observed that Tnt1 actively transposed during somatic embryogenesis, generating an average of 6.37 insertions per line in a population of 19 independent R1 regenerant plants analyzed. In seed-derived progeny of R1 plants, Tnt1 segregated in a Mendelian ratio of 3:1 and no new Tnt1 transposition was observed. A total of 126 flanking sequence tags (FSTs) were recovered from the analyzed R0 and R1 lines. Analysis of the FSTs showed a uniform pattern of insertion in all the chromosomes (1-5) without any preference for a particular chromosome region. Considering the average length of a gene transcript to be 3.37 kb, we estimated that 29 613 lines are required to achieve a 90% possibility of tagging a given gene in the B. distachyon genome using the Tnt1-based mutagenesis approach. Our results show the possibility of using Tnt1 to achieve near-saturation mutagenesis in B. distachyon, which will aid in functional genomics studies of other C3 grasses.

Brachypodium distachyon triphosphate tunnel metalloenzyme 3 is both a triphosphatase and an adenylyl cyclase upregulated by mechanical wounding.

Proteins with a CyaB, thiamine triphosphatase domain (CYTH domain) may play a central role at the interface between nucleotide and polyphosphate metabolism. One of the plant CYTH domain-containing proteins from Brachypodium distachyon, BdTTM3, is annotated in NCBI databases as an 'adenylyl cyclase (AC)' or a 'triphosphate tunnel metalloenzyme'. The divergent nomenclature and the search for plant ACs induced us to experimentally confirm the enzymatic activity of BdTTM3. Based on in vitro analysis, we have shown that the recombinant form of BdTTM3 is a protein with high triphosphatase activity (binding both tripolyphosphate and ATP) and low AC activity. Furthermore, the analysis of BdTTM3 transcriptional activity indicates its involvement in the mechanism underlying responses to wounding stress in B. distachyon leaves.